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Abstract. - The problem of quantum criticality in the context of itinerant ferro- or metam- 
agnetism has received considerable attention [S. A. Grigera et. at, Science 294, 329 (2001); C. 
Pfleiderer et. at, Nature, 414, 427 (2001)]. It has been proposed that a new kind of quantum 
criticality is realised in materials such as MnSi or Sr3Ru207. We show based on a mean-field 
theory that the low-temperature behaviour of the n-layer ruthenates Srn+iRunOsn+i can be 
understood as a result of a Van Hove singularity (VHS). We consider a single band whose Fermi 
energy, Ef, is close to the VHS and deduce a complex phase diagram for the magnetism as a 
function of temperature, magnetic field and Ep- The location of Ep with respect to the VHS 
depends on the number of layers or can be tuned by pressure. We find that the ferromagnetic 
quantum phase transition in this case is not of second but of first order, with a metamagnetic 
quantum critical endpoint at high magnetic field. Despite its simplicity this model describes 
well the properties of the uniform magnetism in the single, double and triple layer ruthenates. 
We would like to emphasise that the origin of this behaviour lies in the band structure. 



Introduction. - The issue of metamagnetism in itinerant electron systems was stud- 
ied both theoretically [1,2] and experimentally [3] long ago. Recently, the investigation of 
metamagnetism in two compounds, MnSi [4] and Sr3Ru207 [5,6], has revived interest in this 
phenomenon, considering it from a new point of view. It has been suggested that these sys- 
tems might display a new type of quantum criticality, connected with a so-called quantum 
critical end point (QCEP), in the vicinity of which the Landau Fermi-liquid theory of metals 
breaks down. In fact the field H„i{T) initiating the metamagnetic transition, the abrupt in- 
crease of the magnetisation, defines a line of first-order transitions in the field (H)-temperature 
(T) plane without any symmetry-breaking. The first-order line ends in a critical end point 
(iJc, Tc). A QCEP occurs if Tc is suppressed to zero as a function of an additional parameter, 
such as the pressure or the chemical composition of the material. 

While theoretical results concerning the properties of the QCEP have been obtained re- 
cently on the basis of phenomenological, low-energy field theories [7], to our knowledge no 
discussion of the microscopic origin of such a QCEP has yet been provided. In the following, 
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we present a model which gives rise to a QCEP of the above kind, which we will discuss 
within a simple Hartree-Fock theory. This model is based on a system with a single, two- 
dimensional electron band where the Fermi level is close to a Van Hove singularity (VHS). 
Our results for the magnetisation as a function of field, temperature and the band filling 
are in qualitative agreement with observed properties of the two- and three-layer ruthenate 
compounds [8,9], and suggest that the general magnetic phase diagram of n-layer ruthenates 
might be understood in terms of band-structure properties. 

Mean-field theory of itinerant uniform magnetism. - We consider a single-band model of 
electrons interacting via an on-site Coulomb repulsion U. In the Hartree-Fock (or mean-field) 
approximation, the free energy per unit cell / is given as a function of the particle density n, 
the uniform magnetisation m (in units of g^B per unit cell) and the temperature by 

fHF{n.,in,T) ^ fa{n,m,T) + Un^ni, (1) 

where = n/2 ± m are the densities of up- and down-spin electrons, respectively and fa is 
the free energy in the absence of interactions, i.e. 

fo{n,m,T)= ^ [na{fia,T) + fia-na], (2) 

where ilo (pier , T) is the grand-canonical potential and /io- (rio- , T) are the chemical potentials 
for the two species of electrons in the absence of interactions. The only information about the 
band structure which enters the mean-field theory is the density of states (DOS) as a function 
of energy, p{lu) ~ <^(<^~£fc)i where Sk is the electron dispersion and the sum is over the 

first Brillouin zone. Given p{u>), one obtains flo{fj,cr,T) = —ksT Jdep(e)ln(l -|- e"''*-^"'^"^) , 
where kg is the Boltzmann constant and /3 = {ksT)"^. The chemical potentials ii„{n„,T) 
are given implicitly by 

■n-a = -df^^^Q = Jdep{e)fT{e - Pa), (3) 

where fxis) = l/(exp(/3e) -I- 1) is the Fermi function. 

In the following, we omit the suffix "HF" in /he- Once the free energy is determined, we 
obtain the thermodynamic equation of state 

h = d,nf = P] - m, (4) 

where h is the external magnetic field multiplied by gps- Equation |0J may be solved for 
m to obtain the magnetisation as a function of n, T and h. The thermodynamically stable 
solution minimises the Gibbs free energy, g — f ~ hm. 

A first-order transition is characterised by a discontinuous jump of the magnetisation, 
whereas a smooth metamagnetic transition occurs if the differential susceptibility x — dm/dh 
has a (more or less pronounced) maximum. The susceptibility is obtained from 

- = C/-EtT^~2C/, (5) 

where ^„(m,T) ^ Jde p^''+^\e)fT{e - p). Note that limT^o An{p,T) = p''"\p), where 
/9'^"^(/i) denotes the n-th derivative of the DOS at the Fermi level. Because dmX\m=o = 
and limm^oo X = 0, a sufficient (but not necessary) condition for either smooth or first- 
order metamagnetism is given by d^x\m=o > 0; which is equivalent to d^f\jn=o < or 
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AqA2 > 3(j4i)^, where An is evaluated at ^ — ^-^ — /^|, i.e. in the absence of magnetisation. 
At zero temperature, this is equivalent to the condition 

PP" > HpT, (6) 

which has been discussed by Wohlfarth and Rhodes [1]. The condition is clearly satisfied if 
the DOS is large and has a strong positive curvature at the Fermi level. 

Such a positive curvature also gives rise to two unusual finite-temperature properties. First, 
x{T)\m=o is not monotonous, but has a maximum as a function of T. Second, the entropy 
s = —drf at low temperatures does not decrease monotonously with the magnetisation, but 
has a maximum at a finite to. Both phenomena are related to the fact that, if p{e) is positively 
curved, the average DOS in a finite region around e^? is higher than the DOS at the Fermi level. 
A sufficient condition for the two properties mentioned above is given by d^d^f\m=T=o < 
or pp" > (p')^, which is a weaker condition than eq. 

Equation © is satisfied if, for example, the Fermi level occurs at a local minimum of 
the DOS. This case probably applies to MnSi [10]. Another example, which wc wish to 
address here, occurs if the Fermi level is close to a VHS. This is the case for the Sr-ruthenates 
which as layered perovskite systems represent quasi-two-dimensional electron systems. In 
two dimensions the VHS is logarithmic, i.e. the DOS has the asymptotic behaviour p{e) = 

In |-^|, where Wi^2 are two parameters with units of energy, which are both of the order 
of the bandwidth(^). For simplicity, we choose Wi = W2 — W in the following, and we extend 
the logarithmic form of the DOS from e = —W to +W. 

One way to study phase transitions is to develop / in powers of to. This has been per- 
formed previously for related systems [1,2, 10], but the application of such a procedure can be 
complicated. At least at zero T, this approach is not successful in the present case, because 
as soon as eq. © is satisfied, 9^/|m=T=o < for all n > 4. In fact, f{m)\T=o is not analytic 
at the specific value of to where one of the Fermi levels crosses the VHS. 

Results. - Despite its simplicity, our model yields a complex mean-field phase diagram. 
In fig. ^ we show the phase diagram as a function of temperature, magnetic field and the 
parameter a;, which is the difference between the actual electron density and the density at 
Van Hove filling, i.e. x = n — 2. 

If the system is very close to the VHS, we find the usual Stoner instability, i.e. a second- 
order transition as a function of temperature. The line of second-order transitions within the 
(r, a;)-plane of the phase diagram is characterised by a diverging susceptibility, or equivalently 
UAo{fJ,{x,T),T) = 1. If the system is pushed further away from the VHS, the simple picture 
of a ferromagnetic low-temperature and a paramagnetic high-temperature phase is no longer 
valid, and we find a weak reentrant behaviour, i.e. the system is ferromagnetic for intermediate 
temperatures, Ti < T < T2. The reason for this behaviour lies in the entropy anomaly 
discussed above. The reentrant behaviour occurs only in a very narrow range of x (0.089 < 
X < 0.091 for U = 0.21^). 

At low temperatures, the (reentrant) transition between the ferro- and the paramagnetic 
phase is of first order, i.e. it occurs due to the crossing of two minima of g{m) = f{m) — hm. 
For example, the quantum phase transition as a function of a; at T = /i = occurs at 
X ~ 0.089, whereas the Stoner criterion Up{£f) = 1 would give a slightly smaller value, 
X ~ 0.081. A metastable ferromagnetic state, which is present in the paramagnetic region close 
to the first-order transition, gives rise to a first-order transition in an external magnetic field. 

(^)As an example, the 7-band of Sr2Ru04 is well approximated by a two-dimensional tight- binding model 
with nearest-neighbour hopping t = 0.17 eV and diagonal next-nearest-neighbour hopping t' = 0.08 eV. For 
this model we find Wi = 2iT^Vt'^ - 4t'2 ^ 1.14 eV and W2 ~ 0.92 eV. 
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Fig. 1 - The phase diagram for U = 0.2W as a function of temperature, magnetic field and x, the sep- 
aration of the electron density from the Van Hove filling. The transition from ferro- to paramagnetic 
behaviour changes from second to first order. The grey surface indicates a first-order metamagnetic 
transition (i.e. a jump in the magnetisation) and terminates at a line of critical endpoints. This line 
is suppressed rapidly to low temperatures and finally reaches zero temperature at a QCEP. The inset 
shows a detailed section of the phase diagram projected on the (a;, T)-plane. 

Fig. 2 - Magnetisation as a function of external field for three different fillings and for temperatures 
fcsT^O, 1, 2, 3, 4, 5- IQ-^M^. 



This transition is is characterised by a jump of the magnetisation without any spontaneous 
symmetry-breaking. The critical field of the transition decreases slightly with temperature at 
fixed values of x. 

The line of first-order transitions (at fixed x) terminates at a critical endpoint {hc{x), Tc{x)), 
which forms a Hne of critical endpoints in the parameter space including x. This line of critical 
endpoints is determined by the equations: d^f = d'^f = and eq. I@J, which must be solved 
for h, T and m at a given value of x. The second derivative is given in eq. (|3J) and an explicit 
expression for the third derivative is = Ai{fii)/AQ{^i) — Ai(^|)/Ao(/i|). 

With increasing separation from the Van Hove filling, the critical point moves rapidly 
to lower temperatures and higher fields. Finally, the line of critical endpoints reaches zero 
temperature, thus giving rise to a QCEP at the point {xqc,hqc)- The low-temperature limit 
is more involved. In fact, the equation d^f — must be replaced at zero temperature by the 
requirement that one of the two spin bands, for example the spin-down band, is at Van Hove 
filling, i.e. /^xIqcep = 0, which implies that Xqc — ^ruqc- As a consequence 9^^/ is not zero 
at the QCEP and actually has a pole. We find 
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Fig. 3 - Zero-field spin susceptibility as a function of temperature for three different fillings, x — 0.095, 
0.1 and 0.11 (from top to bottom). 

2V = {W-2U)e-^, (8) 

where W is the energy parameter in the model DOS. For U — 0.2 W, this yields Xqc ~ 0.29 
and hqc ~ 25 • lO-'^iy. The basic structure of the phase diagram of fig. ^is invariant inside 
the range < U < W/2, although the details depend strongly on U. For U > W/2, both first 
order metamagnetism and quantum criticality disappear. 

Within our model and choice of parameters, the QCEP lies outside the experimentally 
accessible region. Accessing this point would require a magnetic field in the order of 200 
Tesla, and even if such enormous fields could be created, the metamagnetic behaviour {i.e. 
the anomaly in m{h)) is negligible. We find that is suppressed exponentially as the QCEP 
is approached, i.e. Tc{x) ~ exp(a/(a; — Xqc)) with a — (|^)^exp(— w 1.28. This leads 
to critical endpoints at extremely low temperatures for realistic values of the magnetic field, 
which are practically indistinguishable from the true QCEP. The Fermi level, at the critical 
end point vanishes even faster than Tc for x — > Xqc- 

In fig. 13 we show the equilibrium magnetisation as a function of the magnetic field for 
three different fillings and for various temperatures. At x = 0.09, the system is in the reen- 
trant ferromagnetic regime. At very low temperatures, a first-order magnetic transition can be 
found at a small but finite field. A metastable magnetic phase also exists in the absence of the 
external field. True equilibrium ferromagnetism is established only at intermediate tempera- 
tures, and there is a second-order phase transition towards a high-temperature paramagnetic 
state. 

For X = 0.1, the system is in the metamagnetic regime. At low temperatures there is 
a first-order metamagnetic transition. Both the critical field of this transition and the size 
of the magnetisation jump are decreasing functions of the temperature. The discontinuity is 
transformed into a smooth metamagnetic crossover above the critical temperature, fesTc ~ 
1.2-10^'^M^, and at still higher temperatures, the normal paramagnetic behaviour is recovered. 
Further away from the Van Hove filling, at x = 0.11 (fig. He)), the behaviour is similar, but 
the metamagnetic transition occurs at higher fields and Tc is reduced. 

From figs.Elb) and c), we observe that the linear response to low fields is not a monotonous 
function of temperature, but has a maximum at a finite temperature T* which is higher than 
Tc. For fillings close to the reentrant ferromagnetic regime, ksT* approaches 3.2 ■ lO-^W^, 
the temperature where ferromagnetism first appears, as x is lowered. At this point, which 
corresponds to a second-order ferromagnetic transition, the susceptibility is in fact a diverging 
quantity. For larger values of x, the peak in xC^) is much less pronounced and T* increases 
(fig- 13 • A very broad maximum remains even for fillings much larger than Xqc- 
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Comparison with experiment. - In view of the metamagnetic behaviour of the muhi-layer 
strontium-ruthenates (the layered perovskite compounds Sri+nRunOsn+i), it is natural to 
concentrate on the electronic band which is related to the 7-band of the single-layer compound. 
This band, which originates from the Adxy orbitals of Ru, lies close to a VHS. The double-layer 
compound Sr3Ru207 has a modified band structure which preserves two bands with the basic 
7-band structure. Band structure calculations suggest that one band is moved closer to the 
VHS [11, 12]. It is natural to assume that this trend of approaching the VHS in one band 
progresses with increasing number n of layers in each unit cell, so that n is a tuning parameter 
for the variable x in one of the bands. Other aspects, such as a possible change in the ratio 
between the bandwidth and the Coulomb repulsion or more complex multi-band physics, 
are neglected. In fact a phase diagram very similar to fig. ^ where the filling parameter a; 
is replaced by the pressure, was proposed on purely phenomenological grounds in relation to 
Sr3Ru207 [6]. For a given compound, x may also be tuned by external pressure or by chemical 
pressure, such as the gradual replacement of Sr by the smaller isoelectronic ion Ca. 

The magnetisation curves in fig. |5]b) are in good qualitative agreement with the observed 
behaviour in Sr4Ru30io [9] for magnetic fields parallel to the a6-plane(^). Furthermore, fig. HI 
c) is consistent with the experimental findings in Sr3Ru207 [8], where smooth metamagnetism 
has been observed at fields which are 2-3 times larger than the critical field in Sr4Ru3 0io [9]. 
We would like to emphasise that the basic, qualitative properties of the double- and triple-layer 
compounds are in good agreement with our model calculation. Our choice of parameters even 
yields values of temperature and magnetic field which are of the same order of magnitude as 
those found in experiment. From our calculations, we expect that a first-order metamagnetic 
transition should be found at sufficiently low temperatures in Sr3Ru207, but given the rapid 
decrease of as a function of x, the required temperature is likely to be much smaller than 
that attained in the experiment. 

We stress that in our model the QCEP occurs when one band is at the Van Hove fill- 
ing. This circumstance may influence unconventional transport properties in the vicinity of 
the QCEP [5], for non- Fermi- liquid behaviour has been suggested as one consequence of an 
itinerant electron band being at the Van Hove filling [13-15]. 

Instead of increasing the number of layers, an alternative approach to the VHS may be 
through electron doping of the single-layer compound. This was attempted recently by partial 
replacement of Sr by La [16]. The experimental result indeed shows an enhanced magnetic 
response as a function of doping, but neither ferromagnetism nor the typical temperature 
dependence of the susceptibility, with a maximum at finite temperature, was observed. The 
question of the extend to which the doping may lead to an increase in the itinerant electron 
density remains to be answered. If the additional electrons remain essentially localised around 
the dopants, this may yield a similarly enhanced, almost Curie-like, susceptibility [16]. 

Conclusion. - In conclusion, we have presented a mean-field theory of uniform magnetism 
for a single band of itinerant electrons which is close to a VHS. The resulting phase diagram 
as shown in fig. ^ features first- and second-order ferromagnetic transitions, as well as a line 
of (metamagnetic) critical endpoints, which is pushed to zero temperature at a QCEP. Such 
a phase diagram has been proposed on purely phenomenological grounds in the context of 
MnSi [4] and multi-layer ruthenates [6], but no microscopic description has yet been provided. 
Despite its simplicity our model calculation gives a good qualitative description of the magnetic 
properties of two- and three-layer ruthenates [8,9]. The general magnetic behaviour of these 

(^)Sr4Ru30io is actually ferromagnetic below 105K, with an easy axis in the c-direction. Metamagnetic 
behaviour was observed by applying a field parallel to the afe-plane. The very different magnetic behaviour 
parallel and perpendicular to the c-axis is an effect of spin-orbit coupling, which is not included in our model. 
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materials may thus be explained in terms of band-structure properties. A similar scenario 
applies to MnSi, where the Fermi energy is likely to be associated with a dip in the electron 
DOS, and to FeSii-xGcx, which exhibits a transition from a small-gap semiconductor to a 
metallic ferromagnet with increasing x [17]. 

Our analysis reveals the remarkable aspect that the QCEP is characterised by the situation 
in which the chemical potential of one of the bands (minority- or majority-spin) is exactly at 
the Van Hove filling. This should have non-trivial consequences for the metallic properties close 
to the QCEP. In particular, transport properties may show non-Fermi-liquid behaviour [13,15], 
similar to the experimental observation [5]. 

Finally, we comment that, although the model we have presented captures many of the 
qualitiative features observed in multi-layer ruthenates, it remains a relatively crude sim- 
plification. Here we have neglected several aspects, including multi-band effects, spin-orbit 
coupling and disorder effects, which undoubtedly have a certain role in real materials. One 
of the most important questions concerns the effects of spin fluctuations, which are neglected 
within the Hartree-Fock approximation. These aspects are the subject of current studies. 
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